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The form factors near zero momentum transfer are precisely calculated by using the all-to-all 
propagator and twisted boundary conditions. We compare the normalizations and slopes of the 
form factors with ChPT and experiments. 
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1. Numerical simulations 

The matrix element of the K—tnlv decays is parametrized by two form factors 

(K(p')\V,\K(p)) = (p + p'^f+iq^ + ip-p'^f-iq 2 ) (q 2 = {p-p'f). (1.1) 

The normalization of the vector form factor /+(0) is an important quantity for a precise deter- 
mination of a CKM matrix element \V US \ and the search for new physics. For a reliable lattice 
calculation of /+(0), we also examine the consistency of other information in the matrix element, 
namely /-(0) and the form factors' shape, with chiral perturbation theory (ChPT) and experiments. 

In this article, we report on our calculation of //+ _} (<? 2 ) in Nf = 2+ 1 QCD. Chiral symmetry 
is exactly preserved by using the overlap quark action for a straightforward comparison with ChPT. 
At a single lattice spacing a = 0.1 12(1) fm, we simulate three values of degenerate up and down 
quark masses m U( i = 0.015, 0.035 and 0.050 that cover a range of the pion mass 290-540 MeV. 
The strange quark mass is fixed to a single value m s = 0.080, which is very close to its physical 
value m^phys =0.081. We choose a lattice size, (L/a) 3 x (T/a) = 16 3 x48 or 24 3 x 48, depending 
on ntud in order to satisfy a condition M K L > 4 to control finite volume effects. The statistics are 
2,500 HMC trajectories at each combination of m u d and m s . 

We calculate two- and three-point functions 

C p (A*,p) = -4££<^(xV + A0^(x,0), (1.2) 

L 1 x,t x' 

c£ G (A^';p,p') = £ (0Q(x",t + At + At')V^x',t + At)4(x,t)} (1.3) 

L 1 x,t x ",x' 

using an exponential smearing function 0(r) =e~°- 4 l r l for the interpolating operator &J,(x,t) = 
Lr0( r )<?( x + r it)Ysq'{x,t) (P = K orK). We refer readers to Refs. [jl|, |||] for details on how to 
construct these correlators using the all-to-all propagator [^]. 

In order to explore the most important kinematical region q 2 ~0, the meson momentum pW is 
induced by employing the twisted boundary conditions (TBCs) [Q] 

q(x + Lk,t) = e ie q(x,t), q(x + Lk,t) = e~ w q(x,t) (k = 1,2,3), (1.4) 

where k is a unit vector in the k direction. We use a common twist angle 6 in all the spatial 
directions for simplicity. Here we consider the K + — > 7T°/v channel, and impose the TBCs for the 
up and strange quarks. The periodic boundary condition is used for the spectator down quark. Our 
simulation parameters are summarized in Table [I]. 

Table 1: Simulation parameters. We denote the bare quark masses in lattice units by m^ ud s y 



m ud 


m s 


lattice 


M K L 


d 


0.050 


0.080 


16 3 x 48 


4.9 


0.00, 0.40, 0.96, 1.60 


0.035 


0.080 


16 3 x48 


4.1 


0.00, 0.60, 1.28, 1.76 


0.015 


0.080 


24 3 x 48 


4.2 


0.00, 1.68, 2.64 
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Figure 1: Left panel: effective value /o (<7max ; Af , Af ' ) = 2 M K M % R (At, At')/ (M K + M n ) obtained with dif- 
ferent values of At + At' . Data are shifted in the horizontal direction so that the meson source and sink 
operators are located at T/4 — (At + At')/2 and T/4 + (At + At') /2, respectively. Right two panels: three- 
point function Cf 71 for each jackknife sample. We plot data normalized by the Monte Carlo average for jj. =4 
and 1 in each panel. Open and filled symbols are data with and without averaging over the locations of the 
meson source, respectively. 



2. Form factors at simulation points 

We calculate the scalar form factor fo(q 2 )= f+(q 2 ) +/-(<7 2 )<? 2 /C w i at <7max = {Mr - 

M n ) 2 from the following double ratio []|] 

R(At , A _ C 4 ^(A^A^0,0)Cf (Af,Af';0,0) , (M K + M n ) 2 2 2 

{ ' ' Cf K (At ,At',0, 0)C%* (At , At '; 0, 0) a,,a^ J 4M K M n j0[q ™*> ' 

The form factors /{ +j o}(<7 2 ) at q < <7 max are calculated from |6|] 
« Cf n (At,At';p,p')C K (At,0)C*(At',0) J Er±K , g^L n 2,\ UW) n ~ 



Cf*(A^A^0,0)C*(Af,p)C*(Af',p') lM^ + M^ M^ + M^ w; J/ (* 



i max 7 



= (A^A^p,pQC;f*(Af,Af';p,p') _lp 1 E K + E' K 

k C^{At,At';p,P')C^{At,At';p,p') ^ (p+p% (E K -E' x )Ztf)-E K -E' x > { ' 

where Ep 1 (P = 71 or K) represents the energy of the meson P with the momentum pW. Note that 
we can convert / + (q 2 ) to fo(q 2 ) (and vice versa) using the the ratio £ (g 2 ) = /_ ) //+ (<? ), except 
at q 2 msx where R and have no sensitivity to E, (q 2 ). 

The use of the all-to-all propagator greatly helps us identify the plateaux of these ratios. In 
contrast to previous studies with At + At' kept fixed, we can take arbitrary combinations of At and 
At'. As shown in Fig. [TJ, the effective values of /o(<7 ma x) with different values of At + At' exhibit 
good consistency, which gives us confidence about our determination with sub-percent accuracy. 

Another important advantage with the all-to-all propagator is that we can remarkably improve 
the statistical accuracy of meson correlators (and hence their ratios) by averaging over the locations 
of the meson source (x,r) in Eqs. (1.2) and (1.3). As shown in the right two panels of Fig. Ill the 



statistical fluctuation of C^ 71 (pL =4, 1) is reduced by about a factor of 3. This averaging enables us 
to achieve < 1% accuracy for fr + .o}(<? 2 ), and 10-30% for %(q 2 ) 
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Figure 2: Scalar form factor fo{q 2 ) as a function of q 2 . Left and right panels show data at our heaviest and 
lightest pion masses, respectively. We also plot interpolations to q 2 = with various parametrization forms 
together with /o(0) obtained using the pole ansatz (diamond). 



3. q 2 dependence 

We plot our results for fo{q 2 ) as a function of q 2 in Fig. |2[ In this study, we simulate small 
values of |<7 2 | by using TBCs to precisely determine / + (0)(=/o(0)). The data show small curvature 
in our region of q 1 and are well described by any of the following parametrization forms 



/ (<7 2 )=/o(0)(l+C(K7 2 ), Mq 2 )=M0)(\+c q 2 + c iq 4 ), f (q 2 



M0) 



l-ql/Ml 



(3.1) 



pole 



which have also been used in the analyses of experimental data. In this preliminary report, we 
determine the normalization /o(0) and its slope /q(0) using the pole form. The uncertainty due 
to this parametrization is estimated by the largest deviation among the results of the above three 
interpolations. We note that this uncertainty is similar to or smaller than the statistical error. 

The situation is similar for /+ (q 2 ), which is plotted in the left panel of Fig. ||. We observe that 
/+(0) is in good agreement with /n(0) as expected, while the latter has slightly smaller uncertainty 
due to a better control of its interpolation to q 2 = with the accurate data of /o(<?max)- We therefore 
use fo (0) as the normalization of the vector form factor in the following analysis. 
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Figure 3: Interpolation of f+(q 2 ) (left panel) and E, (q 2 ) (right panel) as a function of q 2 
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Figure 4: Chiral extrapolations of /+(0). The left and right panels show extrapolations based on the chiral 
expansion in terms of M^ n k ^/Fq and M 2 K K y/F%, respectively. In the right panel, we also plot /+(0) from 
recent calculations in N f = 2 + 1 [§] and Nf=2 fj] QCD. 



We plot %(q 2 ) as a function of q 2 in the right panel of Fig. |[ Our data show a very mild 
dependence on q 2 with our statistical accuracy of < 30 %. Note also that the leading order (LO) 
analytic terms in the chiral expansion of ^(q 2 ) are independent of q 2 , since ^(q 2 ) vanishes as 
«m s - m u d in the SU (3) symmetric limit [Q]. We interpolate t, (q 2 ) to q 2 = using a linear fit. 



4. Chiral extrapolation 

Figure ^| shows our chiral extrapolations of /+(0) based on SU (3) ChPT. In this preliminary 
report, we employ the fitting form 



/+(0) = l+/ 2 + A/, 
Mj + M 2 I 



h — -z (Hrx + Hkt] ] 



H, 



PQ 



mx 2 F 2 



IMlM 2 
1 + ..,.., In 



M 4 P - 



M 2 



M 2 



■Mi 



F 2 



CQ + Ci 



M 2 +Ml 



F 2 



(4.1) 



(4.2) 



where f% represents the next-to-leading order (NLO) contribution 

(4Mf — M 2 )/3 is used to evaluate Hkt\- We also include the higher order analytic correction Af 



], and the LO relation M 2 



with C{o,i} treated as fit parameters. The Ademollo-Gatto theorem [ JTOp guarantees that fi consists 
only of the chiral logarithms with the single low-energy constant (LEC) Fq, that is the decay con- 



stant in the chiral limit of three flavors. Note that Eqs. ( |4.1| ) and ( |42D can be considered as an 
expansion in terms of M 2 n k ^/Fq. 

We obtain the extrapolation in the left panel of Fig. |] by using Fq = 52.5(5. l) sta t determined 



from our study of the meson decay constants Ql ip. The convergence of the chiral expansion at 
the physical quark mass is /+(0) = 1 — 0.073 [^2] +0.025(6) [Af] in contrast to the conventional 
wisdom that is only few percent correction and Af is even smaller. This is because our estimate 



of Fq is significantly smaller than the phenomenological value Fo = 87.7 MeV[ 12] and enhances the 
chiral corrections ^M 2 ^. K ^/F { 2n . Note, however, that the phenomenological estimate involves large 
N c assumptions, which are not consistent with recent experimental data of the K14 decays [|13|]. 
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The convergence can be improved by switching the expansion parameter to M^jq/^tc as we 
already demonstrated in our study of M n and F % in Nf = 2 QCD [ 14]. The right panel of Fig. |] 
shows the extrapolation using Eqs. ( ]4.1| ) - ( |4.2| ) rewritten in terms of Mj n K \/F%. We obtain a more 
convergent expansion /+(0) = 0.964(6) = 1 — 0.023 — 0.013(6) without a subtle cancellation 
between f% and Af. Note also that recent calculations in Nf = 2 + 1 and Nf = 2 [Q] QCD are 
consistent with this extrapolation. 

Our chiral extrapolation of the normalized slope {r 2 )y n = 6/+(0)//+(0) is shown in the left 
panel of Fig. |5| We employ the NLO expression [0] 



h(x) 



V2L% 
F 2 

1 K 



64n 2 F 2 



-3x 2 -3x+l, . n 1 

— ; In x + - 

2(x-l) 3 U 2 



Ml J \ Ml 

2 



x+ 1 




+ ln 



, (4-3) 
(4.4) 



plus a higher order analytic correction. Note that {r 2 )y n has the NLO analytic term with a LEC 
L r 9 in contrast to /+(0). As shown in Fig. our data are well described by this form and the 



extrapolated value is in good agreement with the experiment [16]. Our estimate L r 9 x 10 3 =4.1(3) 
is slightly smaller than a phenomenological estimate 5.9(4) [O], though the error of our preliminary 
result is statistical only. 

The right panel of Fig. || shows our results for E, (0) as a function of the SU (3) breaking 
parameter M\ —M 2 . In this preliminary analysis, we parametrize the quark mass dependence of 
§ (0) by a simple linear form 



^(0) = do + d^Ml-Ml) 



(4.5) 



which is motivated from the ChPT expression of the leading analytic terms <=c M\ — M 2 Our 
data are well fitted to this form as shown in Fig. |5| We obtain do = —0.006(8) confirming that E, (0) 
vanishes in the SU (3) symmetric limit, as expected. The extrapolation to the physical point yields 
^(0) = -0.094(8) which is consistent with the experimental value -0.125(23) [ph. 



1 See also Ref. jl5| ] for a resummation of effects of sea strange quarks, which could be a source of the large deviation 
between F n and Fq. 
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5. Summary 

In this article, we report on our calculation of the kaon semileptonic form factors. Their 
normalizations and slopes at q 2 = are precisely calculated by using the all-to-all propagator and 
TBCs. We observe a good consistency of {r 2 )^ 71 and E, (0) with experimental results. 

The choice of the expansion parameter is an important issue on the convergence of the chiral 
expansion of /+(0). We note that the small value of Fq also enhances the chiral correction to other 
observables, such as the pion and kaon charge radii [@]. The large deviation F^—Fq can mainly 
come from effects of sea strange quarks, which are difficult to estimate solely from experimental 
data. A more precise determination of Fq is an important task of lattice QCD to assess the reliability 
of the chiral expansion based on SU (3) ChPT. 

Our calculations at different values of m u d and m„ are underway to study systematics of the 
chiral extrapolations. It is also interesting to apply two-loop ChPT formulae to our data. Although 
the formulae are complicated, exact chiral symmetry forbids additional terms due to finite lattice 
spacings and provides us a theoretically clean comparison at the higher order. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 09/10-09). This work is supported in part by the Grants-in-Aid for Scien- 
tific Research (No. 21674002, 21684013, 23105710), the Grant-in-Aid for Scientific Research on 
Innovative Areas (No. 2004: 20105001, 20105002, 20105003, 20105005), and the HPCI Strategic 
Program of the Ministry of Education, Culture, Sports, Science and Technology. 
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